Finite size effects as the explanation of "freezing" in vortex liquids 
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We investigate the effect of thermal fluctuations on the (mean-field) Abrikosov phase. The lower 
critical dimension of the superconducting phase is three, indicating the absence of the Abrikosov 
phase for dimensions d < 3. Within the d = 3 vortex liquid, the phase correlation length £| 
along the magnetic field direction grows exponentially rapidly as the temperature is lowered. For 
a finite bulk system, there is a 3D-2D crossover effect when becomes comparable to the sample 
thickness. Such a crossover effect takes place over a very narrow temperature interval and mimics the 
"first order transition" seen in experiments on clean YBa2Cu307_5 (YBCO) and Bi2Sr2CaCu20g 
(BSCCO) crystals. We calculate the jumps in the entropy, magnetization and specific heat due to 
the crossover and find reasonably good agreement with experiments on both YBCO and BSCCO. 

PACS: 74.20.De, 74.76-w 



I. INTRODUCTION 

Abrikosov's mean-field treatment of conventional Type 
II superconductors in a magnetic field is very accurateo 
because of their extremely narrow critical regime. In this 
approximation, the superconducting triangular vortex 
crystal melts into a resistive vortex liquid via a continu- 
ous phase transition. However because of the strong ther- 
mal fiuctuations about the mean-field solution in high Tc 
superconductors (HTSC) such as YBCO and BSCCO, 
the nature of the transformation between the supercon- 
ducting and the normal phase (the vortex liquid), and 
even the very existence of the mixed phase itself has be- 
came an issue of great interest and complexity. 

Evidence for a first order melting phase transition 
from vortex crystal to vortex liquid has been found 
in recent niagnetization nreasurements on clean crystals 
of YBCOHu and BSCCOcl. The magnetization jump 
AM associated with the melting is found to be about 
3 X lO-^T at 4T for YBCO and 4 x IQ-^T at 5 x lO'^T for 
BSCCO. By measuring the position of the phase bound- 
ary in the H-T plane, and assuming a first order phase 
transition, one can obtain the entropy jump per vortex 
per CuO layer AS" via the Clausius-CIapeyron equation: 



AS _ s^odH,n{T) 



AM 



B 



dT 



(1) 



where B, s and Hm{T) are magnetic induction along the 
c-axis, layer spacing and melting line respectively and 
$0 is the flux quantum. For YBCO, AS is calculated 
to be O.e/cs/layer/vortex at 4Ty. For BSCCO, AS" is 
calculated to be much higher: about 2A:B/layer/vortex 
at 1 x 10~^Tp.and rapidly growing as the temperature 
approaches TJj. Recently, both the magnetization jump 
and the latent heat have been measured on the same 
YBCO crystal. While Schilling et alB confirmed that the 
jumps in entropy/vortex/layer and magnetization satis- 
fied the Clausius-CIapeyron equation, Junod et alB found 



less satisfactory agreement, possibly due to sample inho- 
mogeneities. Both sets of authors report a disappearance 
of the entropy jump in small fields, which again probably 
indicates that the experimental data is being affected by 
sample artifacts. 

Recently reported numerical simulations also favor 
a first order melting transition. Monte Carlo sim- 
ulations using variaj^&i of the 3D XY model give a 
jump in the entropytl'Elil. Sasik and Stroud's 3D Monte 
Carlo simulatipH within the lowest Landau level (LLL) 
approximationllJ (when corrected for an erroneous def- 
inition of AS") and that of Hu and MacDonaldEll yield 
estimates of AS* w O.Gfcs/layer/vortex, in good agree- 
ment with YBCO experiments. However, the results 
from these simulations in the LLL scheme should be 
treated with caution: they were performed using quasi- 
periodic boundary conditions, which imposes an effective 
(spurious) pinning potential on the vortex motioril3. A 
good example of the problem of using the quasi-periodic 
boundary condition is the case of two-dimensional (2D) 
(thin film) simulations within the LLL approximation: 
authors who use quasi-periodic Japupdary conditions see 
first order vortex meltingE3o'Ej'E3EII, whereas simula- 
tions iiuaihich the 2D vortices move on the surface of a 
spheret§E3, which involves no spurious pinning potential, 
see no phase transition at all! There is also no experi- 
mental evidence that thin film superconductors undergo 
a first order melting phase transition. A detailed discus- 
sion of this topic has been given in Ref. 

On the theoretical front, it has been thought for a long 
time that the Abrikosov lattice will melt into a vortex li% 
uid, and that the phase transition will be first orderll3. 
However, there is as yet no detailed melting theory. Over 
the years, many theoretical investigations based on the 
melting scenario have relied upon the Lindemann crite- 
rion: when the spatial fiuctuations of a vortex due to 
thermal excitations become some fraction cl of the vor- 
tex lattice spacing, the lattice melts. Tba Lindemann 
number cl is usually in the range 0.2-0.4c3. It has also 
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been suggested that the mechanism behind the appar- 
ent first order phase trassition is the decouphng of the 
vortex-Hnes to pancakesEJ. However, this does not ex- 
plain the disappearance at the transition of the p,ys- 
taUinc order seen in neutron scattering experimenta23. 

On the other hand, it is well known that the mean-field 
Abrikosov solution is unstable against long wavelength 
thermal excitation of the shear modes of the voctjex lat- 
tice in both the physical dimensions d = 2 and 3E3i2j. By 
calculating the off-diagonal long range order (ODLRO) 
in the low temperature regime, it has been shown that 
the lower critical dimension cIlc of the mixed phase is 
three for extreme Type II superconductors (k >> 1) 
in an external magnetic fielded. (The lower critical di- 
mension of an ordered phase is the spatial dimension 
at and below which a system can no longer sustain the 
long range order associated with that phase at non-zero 
temperature). It was suggested therefore that thermal 
fluctuations will modify the mean-field phase diagram as 
follows: for d — 2, the normal vortex liquid is the only 
thermodynamic phase. For d = 3, there are only the 
Meissner and the normal vortex liquid phases. This the- 
oretical scenario for d = 3 is seemingly at odds with the 
overwhelming experimental and numerical evidence for 
a first order melting transition as outlined above, and 
hence has attracted little attenjtion. 

Recently, however, one of uaHSl has proposed that the 
apparent first order melting transition of a vortex crys- 
tal to a liquid may be a signature of a finite size effect 
rather than a genuine thermodynamic phase transition. 
This idea is based on an extension of Refs. ^ and In 
this picture for d = 3, there exits only one true thermo- 
dynamic phase, the vortex liquid phase, characterized by 
two length scales £\\ and £±; l\\ measures the phase corre- 
lation length along the magnetic field direction, and 
is the range of the (short-range) crystalline order. Both 
^11 and €x are growing exponentially rapidly as the tem- 
perature is lowered, but they only become infinite at zero 
temperature. There is no phase transition to an ordered 
phase at any finite tempera tj.iro . This scenario is that 
of zero temperature scalingt3E3cj. For a finite system, 
the rapid growth of £\\ and £_l as the temperature is low- 
ered has profound consequences. For a bulk sample with 
a slab geometry and the field along the c-axis, the vor- 
tex liquid phase becomes phase correlated along the field 
direction upon cooling when £\\ reaches the sample thick- 
ness Lz ■ Then one has phase correlation right across the 
sample. The behavior of the system then crosses sharply 
over from thai of a 3D vortex liquid to that of a 2D 
vortex liquidEZI. Such a crossover effftct can explain the 
sudden drop in the c-axis resistivitylHa. We will show 
later that this crossover can quantitatively explain the 
apparent jumps in the entropy and magnetization of the 
system observed in experiments. In this simple scenario 
there is no melting phase transition and a vortex crystal 
phase does not occur. 

The rapid growth of the c-axis phase correlation may 
have already been seen in flux transformer experiments 



on clean YBCO crystalsES. Right at the point where the 
apparent first order melting transition occurs, the voltage 
difference between various points on the top and bottom 
of the sample are as if the flux lines moved as rigid rods, 
indicating phase coherence across the sample thickness. 
(We acknowledge that a growing c-axis conductivity can 
also produce the same behavior without any substantial 
degree of phase coherence being preseiit). Rjscent numcp* 
ical simulations in frustrated VillainO, XYQ and LLLE2I 
models have reported rapidly growing phase coherence 
along the c-axis upon cooling. With a conventional first 
order phase transition picture, it is difficult to explain 
the presence of such a growing length scale at the tran- 
sition. In addition, the jumps in the magnetization and 
entropy in Refs. ^ and |3^ are actually rounded as func- 
tions of magnetic field and temperature. The width of 
the transition has been calculated within the crossover 
approached, and is in good agreement with data. If there 
were a genuine first order transition this rounding has to 
be explained on the basis of a sample artifact etc. 



In this paper, we shall follow the zero temperature scal- 
ing idea and focus on the growing length scale and the 
role of finite size effects. The outline of this paper is 
as follows: we start by describing the LJ-Jj approxima- 
tion based on the the work of Eilenbergeio and use it to 
study anharmonic fluctuations in the Abrikosov phase. 
We then reproduce in SectjoiaJllthe result that dLC of 
the mixed phase is threeE3cJcJt3 within the framework 
of the loop expansion and estimate the growth rate of ^|| 
and as a function of temperature in three dimensions. 
We find the 3D-2D crossover line by setting £\\ k, L^, 
and compare it to the experimental melting line with the 
Ginzburg number as the fitting parameter. Both YBCO 
and BSCCO are examined. In Section tV , we argue that 
the three dimensional form of the free energy crosses over 
sharply to the two dimensional thin-film form, mimick- 
ing the sharp changes in the thermodynamic potential 
associated with a first order phase transition. We will 
calculate the jumps in the first derivative of the free en- 
ergy to find the entropy jump/layer/vortex, AScr and 
the magnetization AMcr and compare them with exper- 
iment (where the jumps are usually interpretated as be- 
ing due to a first order phase transition). We will also 
demonstrate that AScr and AMcr satisfy the Clausius- 
Clapeyron equation. With sets of parameters appropri- 
ate for YBCO and BSCCO, we will show that our results 
are in reasonably good agreement with experiments. The 
effect of weak random disorder, which is always present 
even in clean crystals, is investigated in Section ^ Fi- 
nally, we will conclude with a summary and discussion 
Most of the details of the (necessarily) 



in Section VI 



complicated calculations are to be found in Appendices 
A-D. 
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II. THE MODEL 

We start from the Ginzburg-Landau model with a com- 
plex superconducting order parameter \E' for a system 
with spatial dimension d and in an external magnetic 
field Ho along the {d — 2) longitudinal directions. It 
is assumed that the system has an effective anisotropy 
m = (m^, TO^, mil). The free energy functional is 



d^^T 



-ihV - 2e.4i)^'| 



|B-AioHo|' 



2^0 



(2) 



where a cx {T — T^.) and (3 is taken to be a constant. 
The magnetic induction B = V x A is assumed to be 
uniform inside the bulk of the system and parallel to Ho. 
This approximation is valid for an extreme Type II (GL 
parameter k >> 1) superconductor where the fluctua- 
tions in the vector potential A are negligible compared 
to those of the order parameter. Choosing a Landau 
gauge A = B{—y, 0, 0) and restricting the fluctuations of 
the order parameter to the LLL subspace, the free energy 
functional is reduced to: 



2m\\ 



9z 



where /3k = /3(1 — 1/2k^) and an = a + ehfioHo/nn^. 
an = defines the mean-field Hc2 line below which the 
Abrikosov mean-field solution can be written a£3 



*o = aoV'(rlO), 
<p(r|0) = (2ry)i/^expf-^)^?3 



7r(.T + iy) 



C + 



(4) 
(5) 



has zeros which form a triangular lattice with the fun- 
damental unit cell spanned by the vectors rj = (1, 0)io 
and rii = — 1/2,?/ — \/3/2)^o- ^3 is a Jacobian 
theta function and £ — ^Jh/2eB is the magnetic length. 
The spacing of the vortices is given by the flux quan- 
tization condition rji"^ — 2ttI'^ . For a system with vol- 
ume V — LxLyh'^f the number of vortices is = 
LxLy/{2TT£'^). ^'o minimizes the free energy Tlll at the 
value Fmf = -|a_f/pF/2/3„/3A, with ao = Vlaff |//3k/3a, 
and (3a = 1.1596... is the Abrikosov number for the 
triangular lattice. Following EilenbergerEja, we con- 
struct an orthonormal basis set for the A^^-fold degener- 
ate ground states of the operator (— i?iV — 2eA)^, = 
e"''^<p(r|rfe)}, of which ^"0 is a member. Each basis is 
labelled by a vector p = (k, q) where q is the [d ~ 2) di- 
mensional longitudinal vector and k is a two dimensional 
vector confined to the first Brillouin zone (BZ) associated 
with the ideal triangular lattice. The basis states can be 
generated using the relation: Lp{Y\Yk) = e^^''^f{Y + rfc|0) 
for Tfc = (zfe -|- iyk) = (■'^{ky — ik^)- The normalization of 
(p{r\rk) is taken to be |(/?(r|rfe)P = 1, where the overline 
denotes a spatial average over the unit cell. In each of the 
(d — 2) longitudinal dimensions, the allowed values of q 
are integral multiples of 27t/Lz and there are — L^/s 
of them, where Nz is the number of layers. The number 
of allowed rfc is just N^. Therefore, the total number of 



(3) degrees of freedom in a set of {5' p} is N^N-, 



(d-2) 



To sim- 
plify notation, we will drop the bold type on the vectors 
p, k and q henceforth. 

Next we set up the standard formalism for perturba- 
tion expansion around the condensed mod(£3. Since the 
functional in Eq. (j^) is not translationally invariant, the 
thermal average of VE" is spatially inhomogeneous and it 
is convenient to define a spatially averaged quantity a by 
(5^ = I < \1/ > p. At mean-field level, a — ao. Writing 
the fluctuating order parameter around ^'o as 



/(/ci, fc2|fc3, ^4) 



* = Q;(/7(r|0) -I- (^^f 

7f 1.(^1/(0, o|fc,-fc) I 
1 



(5* 



Cpexp(ig.2;).^(r|rfc), 



d r(/j*(r|rfcjv'*(r|rfcj(/j(r|rfc3)v5(r|rfcj. 



(6) 
(7) 

(8) 



cell 



Implicit in the definition Eq. (|8|) is the conservation of 
momentum 5ki+k2,k?.+kii and the transverse integration 
is over the primitive unit cell of area A. /(fci, k2\k^^ k^) 
can be conveniently expressed in terms of lattice sums 
(see Appendix D). Substituting Eq. (0) into Eq. (||), we 
find that the free energy to quadratic order in Cp is given 

by 



where 

= (-\aH\ + 2/3«52/(0, fc|0, k) + , (10) 

Ap^(3^a^\I{0,0\k,~k)\. (11) 
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Eq. (H) can be diagonalized as 

^2 = Tmf + \Y, (i?+|6pp + i?-|ap|2) (12) 



with the eigenvalues E± and their corresponding nor- 
malised eigenvectors e± : 



The soft and hard modes are so-called because of the 
asymptotic behavior for small k: e_(fc) ~ /c* and 
e+(fc) ^ 2 (see Appendix D). The consequences of this 
/c-dependence of e-(fc) for the lower critical dimension 
dhc of the system will be discussed in more detail later. 

We shall introduce a fictitious source field J into the 
LLL functional in Eq. (^ such that 



(e± + 1) - ] 

^ 71 ( -1 



(13) 
(14) 



c.c) . 



(17) 



where e± = 2/(0, fc|0, k) ± |/(0, 0|fc, -fc)| - 1, and ^|| = 
(?i^/2mj||Q!//|)^/^ is the mean-field correlation length in 
each of the {d — 2) dimensions. The form of Eq. ( p^ ) 
implies that the variable Cp can be written as Cp = 
{iap + bp) / \/2 provided a_p = a* and b-p — b* where the 
complex variables Op and bp measure the amount of soft 
or hard mode generated by the thermal fluctuations. The 
number of degree of freedom for each mode is A^^ivi'^"^-'. 
The soft mode and hard mode propagators are obtained 
to lowest order from Eq. (|l^) and are given respectively 

by 



Below Tc, J singles out "^/q oc <i5(r|0) as the condensed 
mode. All results will finally be evaluated at J = 0. We 
rewrite 6"^ in terms of Op and bp as: 



(5* 



where 



bo 

/w 



^(r|0) 



^^(zap + 6p)Qfce"?-Xr|rfc), (18) 



Qk 



/*(o,oifc,-fc) 

2|/(0,0|fc,-fc)| 



(19) 



knT 



Gbip) = < bpbp > = 



knT 



\aH\ie++qq^) 



(15) 
(16) 



The term (fog is real) which renormalizes the amplitude 
of the condensed mode <p(r|0) has been separated out for 
reasons which will be clear when calculating the equation 
of state in Section [II. The sum p implies that except 



(g, fc) = (0, 0), q takes on all momenta in the {d — 2) lon- 
gitudinal space, and k takes on all the permitted values 
in the whole of the first BZ. Expanding Eq. ([l7| ) using 
Eq. (|l^), we obtain the free energy functional 



2Vbo{-&al + a') + ^bli-al + 3a') + ^ + ^ 



/W 8V 



4ei 



aH\ \ a 



2V V 2 



-j(6_(fc) + l)-l+efg2 



2Vboa\Y,W-{k) + l)\ap\ 



-2(e+(A:) + l)-l + dg' 



+ (e+(fc) + l)|6p|^ 



2 r 

-/3«/3a(«\^+6o) E /^(0,pib2,P3)(6;, ^P2^P3 + ^Op^bp^apg bp^Gp^ap^) 
[F{pi,P2\P3,Pi) {a*pyp^ap.,ap^ + ^pi ^p^ &P4 + 4a;^ 5p J 



V 
2V 



-2F{pi,-p2 \p3 , -P4:)bp^ bp^ ap3 a. 



'Pi 



(20) 



where 



F{Pl,P2\P3,P4:) = 5q^+q^^q^+q^Re[Ql^Ql2Qk3QkJ{kl,k2\kz,ki)]. 



(21) 
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In Eq. (|20|), the sum ofpi is over the whole cell subject 
to the constraints on Op and bp. For the practical purpose 
of perturbation calculation that follows, it is most con- 
venient to impose the constraints explicitly and sum pi 
over half over the d-dimensional BZ (i.e. k is restricted to 
half of the two-dimensional BZ and q takes on all allowed 
momenta) . 

The low-temperature free energy functional Eq. (|20| ) 
can be characterised by an effective temperature T — 
/3Kfc_BT/(2f^^^~^|a_ff p), which is related to another pop- 
ular variable ax via f = 47r/|aT|^^^ for d = 3. For 
d = 2, we define fa^ = P^kBT/2eL,\aH\^ = t^I\ocit\^. 
The low and high temperature limits are represented by 
T ^ or aT{pL2T) —oo, and by T ^ oo or aT{a2T 
oo) respectively. Also useful in the following discussion 
is the definition of ar = -{i/Gif^ihty^l^il -h-t), 
where t = T/T^ and h = 5/5^2(0). Gi is the Ginzburg 
number defined in Ref. |2^. Tc and -Bc2(0) are the zero 
field transition temperature and the linear extrapolation 
of Hc2{T) to zero temperature respectively. 



III. THE EQUATION OF STATE 

One of the chief aims of this paper is to carry 
out a loop expansion around the mean-field Abrikosov 
solution — be yoad. the Gaussian approximation previ- 
ously studiedclrc^l — in fact to two loop order. -Ihis loop 
expansion is well known in the 0(n) modeSEj. We 
stress that it is a systematic perturbative approach in- 
volving no ad hoc Ansatz (such as was employed in Ref. 



pi] ). We shall first calculate the equation of state by 
finding what value of a makes < 6o >= 0. Within the 
Gaussian approximation, this just corresponds to putting 
to zero the coefficient of 6o in the functional Eq. ( po|) (see 
diagram Ci of Fig.|l|). This reproduces the mean-field 
results a = ap. 



6o 



FIG. 1. Diagrams for determining the equation of state via 



setting < ho Ci + C2 + = Q up to 0(l3i 
mean-field value of a is given by C\ — 0. 



1/2n 



The 



To one loop order, the equation of state involves the 
tadpole diagrams shown in Fig.|l|. Explicit expressions 
for each diagram are given in Appendix B. The soft and 
hard mode propagators Ga{p) and Gb{p) are labelled by 
a and b respectively. The wavy line in the Feynman di- 
agrams denotes the order parameter 4'o. The thermal 
averaged value of the amplitude of the order parameter 
(5 = (I < ^ > |2)i/2 will be calculated, where the thermal 
and spatial averages are denoted by angular bracket and 
an overline respectively. To one loop order the value of 
a is 



52 



1 - 



-l\aH\{2T:Y 



[ Sk' 


1 d''-^' 


1 BZ J 





-ik') 



-{k')/e 



r,l2 



(22) 



where ^ = Ti /m||. For 3 < d < 4, all integrations are 
finite without cutoffs. At d = 3, the integration involving 
the soft mode (diagram C2) is infra-red divergent. Let 
us concentrate on the singular piece of Eq. (|2^) as d ^ 3. 
Integrating over q' first gives: 



BZ 



d^k' / d''-^' 



1 



e^{k')/q+q'^ 



BZ 
A 



4-d 



.'^"-'dk'. 



(23) 



In the last step, a circular BZ of radius A is as- 
sumed. The integral becomes logarithmically divergent 
at d = 3 indicating that d]^c = 3. Therefore the 
Abrikosov phase is unstable against long-wavelength fiuc- 
tuations in the thermodynamic limit. Such a conclusion 



has also been reached usin g a , similar analysis within 
the harmonic approximationE3c3. The identification of 
e_(fc) « C66^'*fc^/2Q!Q|Q!^f I for small k points to the nature 
of the soft mode which is responsible for the destruction 
of ODLRO at d = 3: it is a long- wavelength elastic shear 
wave. The hard mode, on the other hand, isn associated 
with the compressional mode of the latticecJ. The low 
energy excitations about the ground state, i.e. thcr-spft 
mode, can be described by an effective HamiltonianEll: 



q.keHC 



an 



(24) 



where HC denotes half of the two-dimensional first BZ. 
C66 and ps are the elastic shear modulus and superfluid 
density respectively. The variable Op is just the Fourier 
component of the phase change 9 of the fluctuating order 
parameter 
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1 



(25) 



(Note that the constraint a* = a^p in Eq. ( p5| ) guaran- 
tees that 9 is real.) The soft mode displacements u of 
the vortex lines can be expressed in terms of the, deriva- 
tives of 8: Ux = —(^dBjdy and Uy = (^dBjdxi^. Since 
V.u = 0, the flux line-potions associated with the soft 
mode are shear wavesEj. On introducing the dimension- 
less transverse and longitudinal lengths R — {x/t,y/l) 
and Z — z/^ii respectively, Eq. ( p^ ) becomes in these di- 
mensionless variables 



eff 



ksT 
1 

2f 



,d-2r 



de 



(26) 



where cee = /3KC66/2|a_f/p and ps = m\\PK,Ps/^'^\aH\ 
are the dimensionless shear modulus and superfluid den- 
sity respectively. At the LLL mean-field level, pCcg = 
0.0885 ... and = 1/Pa = 0.8624 . . . respectivelyEa. 

The length scales £\\ and £± over which the ODLRO de- 
cays can be extracted from the singular piece in Eq. (22). 
As mentioned before, we expect that £j_ is also a mea- 
sure of the range of the crystalline order in the transverse 
plane. It will be determined by setting a = on the left 
hand side of Eq. ( p2|) and integrating q e (0, oo) and then 
k G V2{£l\e-^). We obtain 



e±^eexp{A\aT\^^^/2), 



(27) 



with A = 2y'c66Ps = 0.553 . . . Similarly, one can extract 
i\\ by first integrating over k € V2(0,-^~^) and then q 
from 27r(^fr\^7i). This yields 



^11 exp ^^|q!t 



|3/2 



(28) 



As the temperature drops, £± and £|| grow rapidly but 
only diverge in the zero temperature limit {ar — oo) 
where true ODLRO order exists. The functional, forms 
for £± and £|| have been suggested by one of uj£3 using 
a simple renormalization group argument based on the 
effective Hamiltonian Eq. (p6|). However, the number A 
is estimated here for the first timec3. We acknowledge 
that our estimate of A can only be regarded as an or- 
der of magnitude estimate. For example, setting (a/ao)^ 
equal to a constant <1 rather than zero produces a dif- 
ferent value of A. However, we find that the value of A 
quoted here gives an excellent prediction for the position 
of the "melting" line in YBCO (see below). Our calcu- 
lation of the "jumps" in the entropy, magnetization etc. 
is independent of our estimate of A, as we express their 
magnitudes in terms of the measured value of at the 
"melting" line. In Appendix A, we suggest that only a 
non-perturbative approach incorporating the topological 



defects such as entanglements will lead to a quantitative 
estimate of the coefficient A. 

For crystals of the shape normally used in the stud- 
ies of high temperature superconductors, will grow to 
the system dimension before £± reaches the trans- 
verse dimension L^. When this happens, there is phase 
correlation along the c-axis, and the system will then be- 
have as if it were effectively two dimensional. £± is also 
growing exponentially, although slower than £|| , and it is 
expected to be several orders of magnitude times the lat- 
tice spacing at the crossover temperature when « Lz- 
Therefore, when £|| « L^, the system is a vortex liquid 
with quasi-long range order, which explains the appaseujt 
Bragg-like peaks in neutron scattering experimentsr^lri. 
The apparent formation of sharp Bragg-like peaks from 
the rings (expected to be) seen in the structure factor for 
the vortex liquid phase when the temperature is lowered 
is usually attributed to the freezing of the vortex liquid 
to the vortex crvstal phase. However, a recent theoreti- 
cal investigationES has found that such a transformation 
in the structure factor can also take place entirely within 
the vortex liquid phase in the presence of a weak four- fold 
symmetric coupling to the underlying crystal. In fact, the 
angular width of the peaks 59 varies as 69 cx for a 
given coupling to the underlying crystal. This implies 
that the width of the peaks should shrink exponentially 
rapidly when the temperature is lowered. 

The mechanism for the 3D to 2D crossover can be illu- 
minated by a "toy" calculation for the 3D vortex liquid. 
Consider the Hartree-Fock approximation to the prop- 
agator in the vortex liquid phase given in Fig. I This 
calculation is normally done for an infinite system, but 
we shall do for a system of finite width Lz to illustrate 
the 3D-2D crossover mechanism. 



+ 



FIG. 2. The Hartree-Fock approximation. The double and 
single lines denote the renormalized and bare propagators 
G{q) and g{q) respectively. 

The renormalised propagator G{q) (deriving from < 
\I'*(r', z')^'(r, z) >) is related to the bare propagator 
g{q) = ksT / {uh + q^ /2m\j^) such that the "mass" term 

a}i becomes G{q) — kBT/{aji + h^q^ 12171]]) with 



, (3t,kBT ^ 



1 



, P.kBT /2m|| 




(29) 



where the wavevector g„ — 21:111/ Lz for a system of fi- 
nite size Lz with periodic boundary condition. ^|| = 
(?i^/2mj|ajj)^/^ is the renormalized correlation length 
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along the field direction. By using the asymptotic be- 
havior coth(a;) ~ l/x for small x and coth(x) 1 for 
large x, one can see that Eq. (|2^) reduces to the well 
known results for 2D and 3D (Eqs. (23) and (24) in Ref. 
|37| ) in the hmit oi Lz << i\\ and Lz » i\\ respectively. 
Notice that the 2D limit {Lz << of an is dominated 
by just the n = term in the sums. Although ^|| is not 
exponentially growing (because cIlc = 4 in this approx- 
imation), one can see how the behavior of the system is 
controlled by the ratio of the length scale Lz and ^|| , and 
that in the 2D limit one can proceed as if the flux lines 
were straight rods. 

Returning to the full problem, as the temperature 



change required to pass 
^11 « Lz is very smallE3, 



from the regime ^|| <C Lz to 
we believe that the crossover 
has been mistakenly interpretated as a first order melt- 
ing phase transition. Later on, we will calculate the sharp 
step in the magnetization and the entropy/vortex/layer 
due to this crossover and show that it has many features 
of a first order phase transition. 

First, we shall investigate the position in the H — T 
phase diagram where the crossover £|| « Lz takes place. 
In dimensionless units this occurs when ay = where 



-1 2/3 



(30) 



For a typical YBCO crystal of thickness 0.2mm and ^|| « 
lOA, we estimate using our estimated value of A that the 
crossover is at around a^- « —7.9 which agrees witk-,! 
supposed 'melting' line in previous investigationsEj 
For the same sample thickness, but using a typical shorter 
coherence length « 2.0A appropriate to BSCCO, we 
find that ~ —8.9. Note that the position of the 
crossover is only weakly (logarithmically) dependent 
on L^/^ii, whose dependence on T and B will therefore 
be neglected in what follows below. The dependence of 
the parameter on T and B implies that the position 
of the crossover line in the phase diagram should follow 
the power law: 



H 

O 

PQ 3 




82 84 86 88 90 92 

T(K) 



0.04 




70 

T(K) 

FIG. 3. The experimental 'melting' points fitted with 
Eq. (|l|) (solid line) for YBCO and BSCCO. These data are 
read from the bottom panel of Fig. 5 in Ref. ^ and Fig. 4 in 
Ref. ^ respectively. 

Using the YBCO and BSCCO 'melting' data from 
Refs. ^ and ^ respectively, we find reasonably good fits 
with iBo,Ta) « (140T, 92.9K) and {Bo,Tc) w (O.IT, 
94.3K ) respectively (see Fig. ||). By substituting Eq. ^) 
into the definition of a^, we have 



1/3 



BMO} 

Bn 



2/3 



(32) 



Br 



(31) 



where n = 3/2 and Bq is the zero temperature melt- 
ing magnetic field. Strictly speaking Tc is the mean-field 
transition temperature, and fluctuation effects will make 
its value slightly different from the measured zero-field 
transition temperature. 



We can then determine the Ginzburg number Gi of 
YBCO and BSCCO by comparing Eq. (|^and Eq. 
Using i?c2(0) ~ 170T in both casesc3, we obtain 
GiYBCO ~ 0.006 and Giescco ~ 8000. While the for- 
mer is a widely quoted number for YBCO, the later is 
about four order of magnitude bigger than the usual 
value quoted of 0.1. We believe that this might not 
be unphysical for BSCCO, and the argument is as fol- 
lows. What we are using is a phenomenological model — 
the LLL approximation — which is quantitatively useful 
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as long as the effective temperature ax provides a good 
representation of the true temperature and field depen- 
dence. This seems to be the case for YBCO. On the 
other hand, Bo for BSCCO is approximately three orders 
of magnitude smaller than that of YBCO which suggests 
that fluctuations effects are enormous and are likely to 
renormalizc the bare parameters of the theory. Effects 
from the higher Landau level contributions, the quasi-two 
dimensional behavior of BSCCO and the fluctuations of 
the vector potential, which have been neglected in this 
effective model will all act to modify the dependence of 
Q!T on temperature and magnetic field away from the 
simple bare expression for ar as in Eq. (^) . So in prin- 
ciple, one also should expect G«bscco to be a function of 
both temperature and magnetic field rather than a con- 
stant, as we have done above. However, a qualitatively 
useful description of BSCCO may be possible if we are 
prepared to accept a renormalized Gi much larger than 
the bare Gi. Therefore, we will keep an open mind, and 
proceed to investigate what the effective model can of- 
fer in the description of both YBCO and BSCCO. This 
philosophy seems to have been adopted by other authors 
as well. In their recent Monte Carlo simulation, Hu and 
MacDonaldtil had implicitly used a very large Gi in or- 
der to fit their numerical results to BSCCO. On the other 
hand it is possible that the large value of Gi needed to fit 
the data in BSCCO is really telling us that the mecha- 
nism of the two transitions in YBCO and BSCCO are 
quite different e.g. that the crossover idea applies to 
YBCO but that there is a genuine first order transition 
in BSCCO. 

Also of some interest is the angular dependence of the 
crossover line. When the magnetic field is tilted at an 
angle (j> to the c-axis, then the crossover effect will occur 
at ~ Lz/ cosUi)- Using the general scaling approach 
of Blatter et a/.EJ, the dimensionless temperature scales 
as (ignoring as before logarithmic corrections) 

|a^(0)|3/2 w |a^(O)|3/2y^cos2 0-hsin2 4>/T^ (33) 

where F — yyrn\\/rn± is the anisotropy factor. Through- 
out the rest of the paper, we will only examine the case 
of (/) = 0, i.e. a field along the c-axis. 



IV. THERMODYNAMICS 

Although the crossover behavior discussed above is a 
finite size induced effect rather than a phase transition, 
we can still associate with it a "jump" in the entropy 
per vortex per layer AScr and the magnetization AMcr 
due to the narrowness of the crossover region. We be- 
lieve that apparent first order melting signatures like the 
latent heat and magnetization jumps observed in exper- 
iments are in fact due to the entropy and magnetization 
differences between the 2D and 3D vortex liquid. 

Although the jumps in the magnetization and entropy 
are reported to be sharp, they in fact have a finite 



widthfflH. Such rounding of the jumps is usually at- 
tributed to the sample or magnetic field inhomogeneity. 
However, we can explain such rounding as the natural 
width of the crossover effect. A simple prescription to 
estimate the width is to find the small change SaT in 
(set by £|| « L^) required for « 2i^. Using the 
definition of ar, we obtain SaT/oiT ~ {2/3)6B/ B. Sub- 
stituting this into Eq. ( ^0|) and using = 0.2mm and 
^11 w lOA for YBCO, we get SB/B w 0.057 at a magnetic 
field of 4.2T. Using the YBCO magnetization data from 
Fig. 1 in Ref. ||, we estimate that the width is SB « 0.2T 
at B ^ 4.2. This gives SB/B « 0.05, which is good 
agreement with the predicted crossover width. 

If one can calculate the 2D and 3D free energies F2D 
and F^o of the vortex liquid phase, then one would 
naively expect the AS^r = -{s<i>o/BV)d{F3D-F2D)/dT 
and AMcr = -il/V)d{F3D - F2D)/dH. However, 
there are two subtleties involved here and this expec- 
tation is not correct. The first of these is a simplifica- 
tion. From their definitions ar and a2T are such that 

— — \J Lz\a%^rp\^/'^ / Because i^z > C|| for a bulk 
system, a2T is orders of magnitude larger than ar- For 
YBCO, the crossover occurs at d^j, « —7.9. For the 
2D liquid, this corresponds to ~ —1200 (similar esti- 
mates apply to BSCCO). At such a low effective tempera- 
ture, the behavior of the 2D liquid is basically mean-field 
like, and fluctuation effects are negligible. With this in 
mind, the sharp changes in the thermodynamic functions 
between the two regimes can be obtained by subtracting 
the mean-field expression from that of the 3D expression. 
The second subtlety is that not all contributions to the 
entropy or magnetization are sensitive to the effect of 
crossover. For example, the short wavelength contribu- 
tions are not modified when £|| becomes comparable to 
Lz, and so will not contribute to the jumps. Therefore, 
in calculating the 3D entropy, magnetization and specific 
heat jumps, we need to examine all contributions and 
discard the pieces that are continuous over the crossover 
region. 

Before deriving the thermodynamic functions, we 
would like to specify how we envisage infinitesimal 
changes in the magnetic field, i.e. taking the deriva- 
tives of say, the free energy with respect to the mag- 
netic field. We assume a finite system which is allowed 
to change its transverse area L^Ly so that as the mag- 
netic induction changes, the number of vortices inside 
the system, remain constant. The two are related by 
Nfj, — LxLyl^T^i'^ . This framework naturally allows small 
changes in the magnetic induction without the introduc- 
tion of extra vortices, and because of its calculational 
convenience it has also-Jaaea used in the Monte Carlo 
simulations of vorticealJlljS. 

The total 3D free energy of the system can be written 

as 



F = FuF-Nj'4^Q{f) 
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g(f) = -E^^^^f + £^'^'^f^ + (34) 

where Fmf is the mean- field free energy and G{T) is a 
the 3D dimensionless free energy calculated by expanding 



about the mean- field solution, f is a number given by 
the i-th loop contribution to the dimensionless free en- 
ergy (see Appendix C). By definition, the total entropy 
per unit volume is 



S 



1 dF 
VdT 



1 dF, 



MF 



V dT 



N^nkBLz 



AS = 



N^L.kBTw d 



G{f) 
f 



Per vortex per layer, the leading term in the loop expan- 
sion for the crossover entropy jump is 



37r^sfcBg(^) 



1/6 ^*2/3^*-l/3 



4P 



AMrr = 



TTfinHkBtS^^^ df 
2/3A$oe|| dH 



1/6 



GjT) d fT_ 
f Ull 



The first term corresponds to the mean-field entropy 
per unit volume. Across the crossover region, the first 
term inside the square brackets is continuous because 
both Q{T) and d{T/^\\)/dT are continuous. It is the 
first derivative of Q (T) in the second term that gives the 
impression of a discontinuity at the crossover to the 2D 
regime. Therefore the apparent drop in entropy upon 
cooling through the crossover region is the total entropy 
minus all the background pieces (including the mean-field 
contribution) , that are continuous or smoothly varying in 
the crossover region. Thus, the jump in the entropy per 
unit volume due to the crossover effect is: 



(36) 



(37) 



where the superscript * mean that the quantities are eval- 
uated at the crossover. In the last line, we have used 
the approximation 1/t* <C 3/2(1 — /i* — t*) in calculat- 
ing df/dT at the crossover. (This approximation can 
be easily justified by comparing the order of magnitude 
of the two terms using typical YBCO and BSCCO pa- 
rameters at the crossover. In fact one can establish that 
1/r < 3/2(1 -h* - t*) < l/h*). 

In the same way, we can obtain the magnetization of 
the system via the definition AI — ~{l/V)dF/dH . (This 
means that the magnetization is .aot just — < j^'p > as 
in the original work of AbrikosovEJ) . Using the same ar- 
gument as before, the relevant crossover magnetization 
jump arises from the term d[Q{f)/T]/dH. The leading 
crossover magnetization jump is 



(38) 



T^d_ ( g{f) 
i\\ dry f 



(35) 



Again, we have used the approximation 1/ft.*^ 3/2(1 — 
h* — t*) in calculating dT /dH at the crossover. 

It is now easy to see that AScr and AMcr satisfy the 
Clausius-Clapeyron Eq. (|l|) automatically. The gradient 
of the crossover dHcr/dT can determined by differentiat- 
ing Eq. (|3^) directly. This apparent thermodynamic con- 
sistency in the jumps AScr and AMcr has been usedjla 
argue for thej-existence of first order melting in YBCOHu'EI 
and BSCCOo. However, our crossover scenario seems to 
provide a possible alternative explanation. 

How do our expressions for AMcr and AScr compare 
with the actual experimental results? Notice that it is 
the two loop term f in the free energy, rather than 
one loop term i?*^^-* which gives the leading order contri- 
bution to AScr and AMcr- In order to get an estimate of 
AScr and AMcr we need to calculate f^^-*, i.e. evaluate 
the diagrams shown in Fig.0. 






fa £7 £s 

FIG. 4. Free energy diagrams of 0{Pk)- 
Each vertex is O (/?« ), and each contribution from wavy 

— 1/2 

line is proportional to /3k , and so all the above dia- 
grams are of the same order. The single and double ver- 
tex diagrams have an overall negative and positive sign 
respectively, a labels the soft mode propagator, which is 
infra-red divergent at d = 3. The sum of all diagrams in- 
volving the soft mode remains finite although individual 
diagrams are divergent (see Appendix C). The number 
f is estimated to be 4.4 x 10^^. 
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Having found , we can evaluate the orders of mag- 
nitude of AScr and AMcr at some typical field. These 
result should apply to both YBCO and BSCCO provided 
that the appropriate phenomenological parameters are 
used to used to model them realistically. For YBCO, 
we choose s « lOA, ^||(0) « 2.2A and w -7.9. 
This give us AScr- ~ O.Tfcs /layer/vortex and AMcr ~ 
4.0 X 10-^T at 4T. For BSCCO, we use s « llA, 
^11 (0) « I.SA and ~ —8.9, and we estimate that 
AScr « l.Ofcs/layer/vortex and AA-Icr « 0.5 x lO^^T 
at 5 X lO^-iL J'hese results are in good agreement with 
experimeniuoQ. 
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FIG. 5. AScr and AMcr as a function of T for YBCO and 
BSCCO. Data points for AScr{°) and AMcr(*]are read from 
Fig. 4(a) in Ref. ^ and Fig. 5 and 6 in Ref. ^ respectively. 
Solid lines are theoretical fits using Eq. (^ ) and Eq. ( jj^ ) re- 
spectively, with material parameters discussed in the text. 

More importantly, our crossover results predict that 



AScr - h*-'^/^ - (1 - f*)-i/2 and AMcr - h*^/^ - 
(1 — as T approaches Tc. Fig.|| shows the tem- 

perature dependence of AMcr and AScr compared to 
the experiments using material parameters as previously 
mentioned. ThCjYBCO and BSCCD data points are read 
from Wclp et a/.B and Zcldov et alu respectively. Our re- 
sult for AMcr not only agrees reasonably well with the 
bulk of the YBCO and BSCCO data, it also accounts for 
the general temperature dependence quite well. The di- 
vergence of AScr near Tc is consistent with the BSCCO 
data. However, thisjs apparently at odds with the YBCO 
data, as Welp et alB see vanishing AScr as T — » Tc- The 
authors themselves suspect that this is due to the influ- 
ence of sample inhomogeneity very near TcEI. In view 
of this, more weight should perhaps be given to the low 
temperature points when comparing our results with the 
data in Fig. Q We believe that if sample artifacts could 
be removed, the divergence of AScr as T ^ Tc in YBCO 
would be revealed. 

A consistent picture of even the form of the crossover 
line is not available from experiments. Zeldov et al.a 
deduced an exponent n w 1.55 for BSCCO, which is 
close to the results expected from the crossover meclt 
anism (Eq. (31U. On the other hand, Liang et al.a 
and Welp et alB deduced a YBCO melting line with a 
smaller slope (n ~ 1.34, 1.36 respectively). The confu- 
sion is further esemplified by the YBCO measurement of 
Nishizaki et alEJ where a power law appropriate to the 
London model (n « 2) was used. They estimated that 
A^cr ~ &kB at IT and 25kB at 3T. Both the magnitude 
and the temperature dependence of their results are in 
total disagreement with our results here and other YBCO 
measurementgj'ErB. 

In our approach, there is no qualitative difference be- 
tween YBCO and BSCCO. They just differ in having 
vastly different effective values of Gi. But our treatment 
for BSCCO requires the insertion of a value of Gi hugely 
renormalized by fluctuation effects. Gi would then be 
expected to be a function of temperature and field, and 
not a constant as assumed here. Inserting such tempera- 
ture and field dependence (if known!) could improve the 
fits in Fig. H Without inclusion of disorder, our model 
does not explain the vanishing of AScr at a lower tem- 
perature critical point as observed in BSCCO by Zeldov 
et aln (neither do the melting nor decouphng models). 
Such behavior is thought to be disorder induced. The 
effect of random disorder will be discussed in more detail 
in Section^, and seems consistent with the data of Ref. 

i 

The receuLj-. advent of reliable calorimetric 
measurementa23'El has also made available specific heat 
data for comparison. This has motivated us to calculate 
the leading crossover value for the specific heat jump. 
The total specific heat capacity is given by 
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C _dS _ 1 d'^FMF N^L.ksn 
T ~ df~ ~V dT^ ^ 2(3aV 



g{f) fT 



Td^( git)\ d_ fT\ d_ ( g{f) \ 



(39) 



The first term is just the mean-field part of the specific 
heat. The first term inside the square bracket in Eq. ( |39| ) 
is continuous through the crossover region. Only the last 
two terms with the derivatives of Q are sensitive to the 
crossover and hence contribute to the crossover specific 
heat jump ACcr. After some algebra, we find that to 
leading order in the loop expansion, per vortex per layer. 



T 2/3^^11 (0)r,|aT I 



1/2 



(40) 



In order to make it more convenient to compare with 
the experiment of Schilling et a/E3, we convert this re- 
sult to units of mJ/moleK^ using the scale in Fig. 4(b) of 
fnamely, 1 fc^ /vortex/layer = 0.6 mJ/moleTK), 



Ref. 30 



so Eq. ( io|) in the units of mJ /mole T is given by 
97r2f(2)s5^2(0) 



0.6 X 



2/3^^11 (o)rc 



1/2 



(41) 



Using the appropriate parameters for YBCO, we found 
that the leading ACcr/T is constant at 0.3 mJ/moleK^, 
which is approximately four times smaller than the data 
in Fig. 5 of Ref. ^ (see Fig.^. This discrepancy may 
be due to our neglect of the higher order terms in the 
loop expansion. In any case, as a first order approxima- 
tion, our results give a useful qualitative description of 
the current experimental data. 



2.5 



a 1.5 

H 

0.5 





1 1 1 1 




1 

A 








A<g, 


o 


o 






O 
8 A 




A 

O 


- o 


o 






O 

o : 

o 















76 



78 



80 



82 



86 



88 



90 



92 



84 

T(K) 

FIG. 6. Specific lieat diflFerence due to tlie crossover from 
3D to 2D. Data points witli magnetic field parallel to the 
c-axis are read from Fig. 5 of Ref. ^ The circular and trian- 
gular open symbols represent ACcr/T based on specific heat 
and magnetization measurements respectively. The solid hor- 
izontal line is our result using Eq. (klh using a set of parame- 
ters appropriate to YBCO. 



V. EFFECT OF RANDOM DISORDER 



In this section, we investigate the effect of quenched 
disorder (which is always present even in high qual- 
ity YBCO and BSCCO crystals), on the order in the 
Abrikosov phase. The effects of quenched disorder on 
the "melting" transition of YBCO and BSCCO at high 
magnetic fields is well documented. We will just mention 
a few of them here. Safar et alxB reported the existence 
of an "upper critical point" in untwinned YBCO at a high 
magnetic field beyond which the sharp drop in resistiv- 
ity disappeared. In a recent report ofj-Lhe specific heat 
measurements on YBCO, Roulin et alJt3 have suggested 
that the termination of the "melting line" takes place 
at about 14T. Using a sensitive local Hall probe, Zeldov 
et ala also reported similar feature in BSCCO, albeit at a 
much lower magnetic field (0.038T). It is widely believed 
that at high magnetic fields the pinning of the vortices 
by disorder is more cff pc ij i^ fc. and the "first order melting 
transition" is removedcaE^u. In an illuminating experi- 
ment, Fendrich et aLc3 have directly demonstrated that 
the suppression of the sharp kink in the resistivity drop in 
YBCO is a disorder-induced effect. They measured and 
compared the resistivity before and after a controlled in- 
troduction, by electron irradiation, of point defects in 
the sample. Furthermore, they showed that the sharp 
kinks in the resistivity drop can be recovered by reduc- 
ing the density of the ipjut defects through subsequent 
annealing of the samplecS. At a low magnetic field, dis- 
order is also thought to affect the properties of YBCO. 
The existence of such "lower critical point" produced by 
disorder or sample inhomogeneity (by an unknown mech- 
anism) was recently invoked to explain the disappearance 
of the jumps in magnfttization and cntropy/vortex/layer 
in YBCO experimentsoB'El at low fields. 

With this motivation we proceed to investigate the ef- 
fect of disorder on the Abrikosov phase. For simplicity, 
we will neglect the effect of thermal fluctuations and focus 
on very weak spatially varying disorder characterized by 
a locally varying transition temperature T(r). As aJ 
approximation, we adopt the "random Tc" approachE; 
i.e. r(r) is assumed to have a Gaussian distribution: 



< r(r) > = 0, 
W 

<T(r)T(r')> =Y(5(r-r'), 



(42) 



where ^ • • • ^ denotes averaging over all T(r) configu- 
rations. The GL functional in the LLL approximation is 
given by 
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(a// +T(r)) l^-l 







2m|| 


dz 



(43) 



Our plan is to take the disorder to be weak and investi- 
gate the effect of the disorder on the pure system ground 
state vJ/q = aoip{r\0). Expanding in terms of Op and bp 
about this, we obtain the functional up to quadratic or- 
der 



MF 



p 



(44) 



where J^mf is the ground state free energy. The variables 
ap and bp are now disorder, rather than thermal, induced 
fluctuation amplitudes in the soft and hard modes re- 
spectively. At this order, the minimum free energy is 
determined by the conditions 



dXr, 



0, 



(45) 



where Xp £ {Re[ap], Im[ap], Re[bp], Im[bp]}. A measure 
of the effect of the importance of the disordering effect is 



1 



g.feeHC 



> 



q,k£HC 



/(0,fc|0,fc) - \I{0,0\k,-k)\ 
{2\aH\e^+iq^y ' 



(46) 



If S is infinite it means that the disorder has per- 
turbed the system so much that the nature of the low- 
temperature state is completely altered by it. Integrating 
over q € (0, oo) in Eq. (46) first, and concentrating in the 
long- wavelength fiuctuations (small fc), we obtain 



22rf- 



W 



dkk 



2d~9 



(47) 



where W = W/(2\aH\ ^ ) is a convenient measure 
of the strength of the disorder analogous to the dimen- 
sionless temperature T defined earlier. Notice that the 
integral for 5 is infra-red divergent below d = 4 as e — > 0, 
implying that weak disorder will destroy the Abrikosov 
phase at and below four dimensions. One can easily show 
that a similar expression like Eq. ( ^6| ) for the hard mode 
is finite, for d < 4. 

Since the lower critical dimension (Ilc is four in the 
presence of disorder, it is expected that the crystalline 
order in the three dimensional vortex liquid will have a 
power law dependence on the strength of the disorder. 
What then is the range of this crystalline order i?x as 
a function of Wl The ordered phase will only exist if 5 



is small. By setting « 1, and inserting a lower cutoff 
in the integral in Eq. (^) corresponding to the smallest 
wavevector which can be associated with the crystalline 
order in the vortex liquid phase, i.e. e = (27r)^/^£/i?j_, 
gives for (i = 3 



3 -3/2^1/2 \ 1/2 



^66 



W 



(48) 



The range of the c-axis phase correlation R\\ can be 
obtained likewise by evaluating the integral on the right 
hand side of Eq. ( |47[) over k G (0, A) first, and then inte- 
grating over q G 27r(^|7^, i?|r^). This gives 



2'^7rc66 
W 



(49) 



As expected, both R±^ and R\\ are growing alge- 
braically, rather than exponentially, as a function of W ^ 
and are only infinite in the limit of zero W . The proce- 
dure used to obtain R±^ and i?|| arc similar in spirit to 
that of the treatment of weak disorder by Larkino. 

In order to investigate in detail how the disorder mod- 
ifies the crossover effects discussed earlier, one would 
have to take into account both temperature and disor- 
der induced fluctuations in the perturbation expansion 
about the mean-fleld solution, which is a complicated 
task. However, we can get a qualitative idea by sub- 
stituting the crossover values T* and B* into Eq. ( |4^ ) 
and (^^. This would be expected to give, to the lead- 
ing order, the effect of disorder on the crossover. At 
the crossover, we take cx 1/Bcr, ^\\ oc \aH\^^^^ and 

{anl (X {1 — h* — t*) w bII^ . This gives the magnetic 
fleld dependence of the order in the vortex liquid phase 
at the crossover as 



Ri_ cx l/{B,rWf/^, 
Rn CX 1/{WBU'). 



(50) 
(51) 



Therefore the order decreases as one moves along the 
crossover line from low to high magnetic field. Our sce- 
nario of the effect of disorder is that when it is weak, so 
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that i?|| > Lz, it plays little role in the thermodynamic 
properties of the superconductor. If R\\ < L^, it takes 
the role of Lz in our previous crossover calculation. For 
strong disorder, may be so small that the sharpness 
associated with the crossover is removed and the "jumps" 
disappear. This at least seems to explain the existence 
of an "upper critical field" . 

VI. CONCLUSION 

We have shown that within the framework of the loop 
expansion, the Abrikosov phase is destroyed by thermal 
fluctuations at and below three dimensions in the thermo- 
dynamic limit and the only thermodynamic phase above 
Hci(T) is the normal vortex liquid phase. However the 
range of the ODLRO, which is characterised by £± and 
£|| in 3D, is growing exponentially upon cooling and di- 
verges only in the zero temperature limit. We calculated 
the growth of ^|| and £± within the loop expansion and ar- 
gue that the apparent sharp features seen in YBCO and 
BSCCO specific heat and magnetization experiments are 
actually due to the crossing over of the fluctuation behav- 
ior from 3D to 2D when becomes comparable to the 
system thickness. This is in contrast with the widely held 
belief that there exists a genuine first order vortex crystal 
to liquid melting transition well below Hc2{T) line. We 



demonstrated that the entropy/vortex/layer AScr and 
the magnetization jump AMcr due to the crossover sat- 
isfy the Clausius-Clapeyron equation without invoking 
the presence of a first order phase transition. We also 
show that AScr and AMcr can give a reasonable account 
of the magnitude and general temperature dependence of 
the jumps in YBCO and BSCCO. The only free param- 
eter, Gi, is obtained by fitting the experimental "melt- 
ing" line. Our estimate of the jump in the specific heat of 
YBCO is also of the same order of magnitude as the-post 
recent specific heat measurements of Schilling et alEB. Fi- 
nally, we have investigated the effect of quenched short- 
range disorder on the Abrikosov phase. The lower crit- 
ical dimensions d^c is found to be four, implying that 
the short range order in the 3D vortex liquid phase has 
a power law growth as the strength of the disorder is 
reduced. We also demonstrated that random disorder 
tends to remove the sharp crossover effect as the mag- 
netic field is increased, providing a possible explanation 
for the existence of the upper critical field. 
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APPENDIX A: SHEAR MODULUS AND SUPERFLUID DENSITY 



In this section, we calculate the one loop correction to the dimensionless superfiuid density ps and the elastic shear 
modulus C66 ■ We assume that d = 3 + e where e is an arbitrarily small number so that the Abrikosov lattice exists at 
low temperature. The limit e — > will then be taken to get results relevant to the bulk c? = 3 system. Both cgg and 
Ps to one-loop order can be extracted from the equation of state Eq. (j2^) and the renormalized soft mode propagator 
= + + 0{(3l), which involves the two leg diagrams shown in Fig.^. 




Mr 



A^2 ^3 

FIG. 7. Two leg diagrams of order 0{/3k)- 



Each of the above diagrams is of 0{(3k) since each vertex and wavy line contribute 0(13^) and 0(l/\fPZ) respectively. 

^ oS^^f, we have 



The full expression for the M.i and the calculation of ps is given in Appendix B. Writing /Os = 1 / fi"^ — p. 



n (1) 



2<i+lj^rf/24 



d^k 



(d-2)/2 , (d-2)/2 



(d-4)/2 



BZ 



(d-2)/2 (d-2)/2 



(ci-2)|/(0,0|A:,-fc)| 



(Al) 



n (1) 



0.295. 



For d 3, Eq. (Al) can be estimated numerically to give p. 

In principle, the one loop correction to cgg can be extracted from the propagator Eq. (B7) in a similar manner by 
setting q — and expanding in k about k ~ 0. However, it is more convenient to start from the definition of the 
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elastic shear modulus. Following LabuschcS, we allow the distortion of the ideal triar 
that the area of the primitive cell of the first BZ is preserved, cge is then defined as£ 



{ular lattice with the constraint 



C66 



2 d^F{r]) 



V dif 



(A2) 



where is a dimensionless variable specifying the shape of the unit cell with area 2"^^ (see Appendix D) and Firf) 
is the free energy as a function of 77 as in Eq. (p^). However, only the one loop free energy is needed here. The 
second derivative of F is evaluated at the value of i] appr op riate for a triangular lattice, denoted by 'riA{= a/3/2). The 
mean-field shear modulus has been obtained by LabuschE3, who found cgg^'^'' ~ 0.0885. By extending this method to 
one loop order we have: 



C66 



5(0) 



5(1) _ 



(A3) 
(A4) 



where 



(7(77, rf) = / d'k 

JBZ 



^(d-2)/2 ^ ^(d-2)/2 



(A5) 



Our results show that a{d,r]) has a maximum at r]/\ for any d (see Fig.^). This is consistent with the definition 
Eq. (A2) where the second derivative is evaluated at a saddle point, i.e. at rj^. However, the total free energy would 



still have a minimum at rj^ with the curvature of the free energy slightly altered since the mean- field energy F]\{f 
dominates for small T. For d ^ 3, we found using the data in Fig. ^ that cge''^'' ~ —0.0484. 

Because of the destruction of the Abrikosov lattice at d = 3, we might have expected that both cge and jjs should 
be singular. However, at one loop order both of these quantities are found to be finite in the limit d 3. The 
perturbative one loop corrections to ps and cge are not 'aware' of the destruction of the lattice at d = 3 in the 
thermodynamic limit. 

However, we believe that in a nonperturbative approach, which would include topological defects such as entangle- 
ments, both Cge and ps would vanish. Notice that in the (2 -I- e) treatment of the 0(n) modeJ, n > 2, the mechanism 
of the transition is the vanishing of the superfluid density, rather than the order parameterEJ. In our case, neither ps 
nor C66 is driven to zero by small amplitude thermal fluctuations, and the mechanism of the transition is the vanishing 
of ODLRO. 



1.435 



g 

4n' 



1.43 



1.425 



1.42 



1.415 



d = 2.98 J 



ooooqooooo, 



" d = 3.00 
° □ d = 3.02 
. d = 3.04 



d = 3.06 J 



0.8 0.82 0.84 0.86 0.88 0.9 0.92 0.94 

FIG. 8. The function a/in'^ as a function of rj for different d. The vertical dashed line denotes the position of the peaks at 



77A. 
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APPENDIX B: THE SOFT MODE PROPAGATOR AND EQUATION OF STATE 



In t his appendix, we give the exphcit expressions for the one loop diagrams discussed in Appendix A and in 
Section III, Note that the expressions for M^lp) are in abbreviated form, i.e. they are not explicitly even function of 
the external momentum p = (k,q). Such symmetry must be imposed. The primed parameters p' — {k',q') denotes 
the internal momentum that has to be integrated over. We shall also use the notation J dp' = J d!^'^ q' Cpk' / {2ttY . 



C2 

C3 
Mi{p) 

M2{p) 

Msip) 
Mi{p) 
M^ip) 



7 

2f3^(3A&kBTVv 
7 



dp' 
dp' 



e_(fcO + l 
^e^{k')/q + q 

e+{k') + l 
e+ik')/C^ + q 



ksT 

keT 

keT 

keT 
4/3^ lag I /3a 
iksTr 



dp'Gaip') {Hk', k\k', k) + 2Re [QlQfl{k, ^k\k', -fc')] ) 
dp'Gbip') {I{k',k\k',k) - 2Re [Ql,Qfl{k,-k\k',-k')\) 
dp'Gaip') + l + 

dp'Gbip') + l + 

dp'Gaip' ~p)Gb{p') X 

{ |/(0, k'\k, k' - fc)| V |/(0, k\k', k - + |/(0, k - k'\k, -fc')|^ 
+4i?e [QlHO, k\k', k - k')I{0, k - k'\k, -k')] 
-ARe [g^/(0, k'\k' - fc, A:)/(0, k' - fc|fc', -k)] 
-4Re [Ql_k,I{0, k'\k' - k, fc)/(0, k\k - k', k')] 
+8Re [Ql^QlQl^kl'iO, k'\k, k' - k)] 
+8Re [Ql,Ql^Ql_,,I^{0, k\k', k - k')] 
+8Re [QlQlQlL.l'iO, k' - k\k', -k)] 
+8Re [QlQlLkH^, k - k'\k, ~k')I{k ~ k', k'\0, k)] 
~8Re [Q*f3Qll{Q, k'\k' - k, k)I{k - fc', A:'|0, k)] 

-8Re [Ql?Ql,_J{0, k'\k, k' - k)I{k, -fc'|0, k - k')] }. 
The renormalized soft mode propagator to one-loop order is such that 

G-'{p) = G-\p) +Mi{p) + M2{p) + Msip) + MM + M,{p). 



(Bl) 

(B2) 

(B3) 
(B4) 
(B5) 
(B6) 



(B7) 



(B8) 



The propagator Ga{k) is still massless. This can easily showpr,by setting the external mpmcntum p = 0. Indeed, Ga 
is expected to be massless to all orders in the loop expansioncJ. The superfluid densityE3 can be calculated using the 
following relation for small q: 



Substituting Eq. (|2|) into Eq. 



PsQ 
|), we have: 



2kBTGa\q,k = 0)a 



(B9) 



P^PAksT 



dp' X 
e.{k') 



\ e+ik') + 2 A\aH\\I{0,0\k',-k')\ 
\e+ik')/q + q'^ + 7(e+(fc')/eff + <7'2)^ 

32\aH\h'^\I{0,0\k',-k')\' 



+9" lHd-2){e_ik')/q + q'^){e+{k')/q + q'^)3 



(BIO) 
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All integrals in Eq. (BIO) are infrared convergent. On integrating over q' , Eq. (BIO) yields Eq. (Al). Similarly, one 
can calculate cge in principle via 



C66^4^4 ^ 2kBTG-\k -> 0, g = 0)a2. (Bll) 

However, expressing Mi in terms of explicit polynomial of small /c is a challenging task. One would expected that the 
leading order should be fc'* corresponding to a dispersion appropriate to the elastic shear mode. A more straightforward 
way of calculating the correction to cge is to use the definition discussed in Appendix A. 



APPENDIX C: LOOP EXPANSION OF THE FREE ENERGY 



As discussed in Section [V, the starting point of our calculation of the jumps in the magnetization and the entropy 
is the 3D free energy. In this appendix, we derive the 3D free energy expansion about the mean-field solution to two 
loop order starting from the definition: 



F = -/csTln m J [Vap][Vbp] exp(-.F/fcsT) 



(CI) 



At the one loop level, we substitute the Gaussian functional Eq. (|lj) into Eq. (|Cl|) , and we obtain the correction to 
the mean-field as: 



p{i) ^ keTy. 

2 ^ 

k,q 



In 



2|aff|e+(fc)+7g2 



In 



2\aH\e-{k) + 



2 \ 1 



47^3 



27rkBT J V 2TTkBT 

d2fc(e+(fc)('^-2)/2^^_(^)(<i-2)/2) 



k.q 



In 



79 



(C2) 
(C3) 



The first and second term in Eq. (C3) correspond to F'^^^ — F'^^\aH = 0) and F'^^^(aij = 0) respectively. The second 
term is dependent on an ultraviolet cutoff in the longitudinal vector q. The cutoff is of the order of the reciprocal of the 
layer spacing of the model. Since we are interested in the long wavelength limit, we follow the standard prescription 
of absorbing this term into the normal phase energy (see Ref. Using this results in Eq. (^), with S^^'^ « 0.524 
for d = 2>. 

For the sake of simplicity, we set d = 3 at the outset when discussing the two loop calculation. The total two loop 
free energy can be written as 



f(2) 



(C4) 



£s is the sum of energy diagrams involving the soft mode and £h is the sum of diagrams containing the hard mode 
only. After integrating over the longitudinal component q, and expressing the remaining transverse integral in terms 
of the dimensionless BZ (to simplify notation, the tilde on k and k' are dropped), we have 



£3 
£a 



5 

E- 

1=1 



d^kd^k' ^f+{k,k') 



(47r2)2 ^e_(fc)e_(fc') 
d^kd^k' -2/_(fc,fc') 



(47r2)2 ^e^(k)e+ik') 

d^kd^k' (£_(fcO + l)(£_(fc) + l) 

(47r2)2 2Ve_(fc)e_(fc') 
d^kd^k' (e+(fc') + l)(e-(fc) + l) 



(47r2)2 
d^kd^k' 



v/e+(fc')e-(fc) 



gik,k/) 



(47r2)2 ^e+{k + k')e.{k)e.{k'){^e+{k + k') + ^M^) + ' 



(C5) 
(06) 
(07) 
(08) 
(09) 
(CIO) 
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where 



f±{k, k') = I{k, k'\k, k') ± 2Re [Ql?Qll{k, -k\k', -k')] (Cll) 
g{k, k') = ime{Ql-yQk'QkI{Q, k' - k\k', -k))f 

+ 16{i?e[QI.Qfc.gfc_fc./(0, k\k', k^k')+ Ql,QkQk-k'IiO, k'\k, k' - k)]Y 

- 32Re{Ql_k,QkQk'I{0, k-k'\^ k', k)}Re{QkQk'Qk'-kI(0, k\k', k - k')} 

- 32Re{Ql_f,,QkQk'IiO, k-k'\- k' , k)}Re{Qk'QkQk'~kIiO, k'\k, k' - k)}. (C12) 

Each integral involving Xj ^JTZ is logarithmically divergent in the infrared-limit, but the sum of the integrands of 
Ei is finite in the infra-red limit (fc', k — > 0), and hence Es is finite. Numerically, we find that Es « 0.8. The diagrams 
in E\i are given by the expressions: 



Eh 

^6 

^7 

^8 



i=6 



<P-kd^k' -U{k',k) 

1 



-0.89 



fk 
d^kd^k' 



2.07 

h{k,k') 



where 



(4^2)2 + k)e+{k')e+{k)ye+{k + k') + ^M^O + VMfc)] 

/i(fc, fc') = 4i?e[Q^+fc,/(0, fc'l - fc. A: + fc')/(0, /c| - fc', /c fc')] 
+8Re[QfQl,I{0, k\ -k',k + k')I(-k, k' + k\0, k')] 
+16 [ReiQl,QkQk+k'HO, k'\ -k,k + k'))f . 

All integrals can be evaluated numerically, and we found that E^ « 1.80. 



0.62, 



(C13) 
(C14) 

(C15) 
(C16) 



(C17) 



APPENDIX D: THE FIRST BZ AND THE FUNCTION I{Ki,K2 \ Ka,Ki) 

In this appendix, we briefly outline a procedure for evaluating integrals involving /(fci, k2 \ fca, k^) over the first BZ. 
Intgerals over the longitudinal vector q can be done analytically, but the integration over k in the first BZ is done 
numerically. ■— ■ 

Restricting ourselves to the class of centered rectangular latticed^, the fundamental unit cell of the vortex lattice 
can be characterized by two primitive vectors rj = (l,0)^o ^ii = (1/2, '7)^0: where Iq is the spacing between 
vortices. The flux quantization condition determines the area of the unit cell as 2ttP{— rji"^). The corresponding first 
BZ has an area of Ai:"^ /rji"^. It is convenient to rescale the transverse length by \/27r^, and construct a dimensionless 
unit cell with primitive vectors r[ ~ (1/y^, 0) and rjj = (l/2y/^, ,Jf[), giving an area of unity. A dimensionless vector 
k can be defined as \/2/Klk. The corresponding dimensionless BZ in the reciprocal lattice has area 47r^. By changing 
7y, one can construct a unit cell of different shape with the area remaining unchanged. The ideal triangular and square 
lattices correspond to 77A — '\/3/2, an d r\ u ~ 1/2 respectively. In this paper, all the integrals are evaluated at t^a 
except in the calculation of cge in Eq. ( |A2| ). i— ir— 1 

In general, the function 7(^1,^21^3,^4) can be expressed in terms of gauge invariant reciprocal lattice sumsE^L^. 
To simplify notation, we will ignore the tilde assigned to the dimensionless k vector and exploit the conservation of 
momentum ki + k2 = ka + k4 on each quartic vertex. 

/(fci,/c2|fc3,fc4) = exp [-|k2 - k4p/47r- i{ky2 - ky4){kx2 - kx3)/2Tr] 

6m,n exp [-{ky2 + ikx2){igx - gy)] 

X exp [-{gyky4 + g^k^^) + i{gxkyz - gykxz)] , (Dl) 
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where [gx , gy) = {m^, (n — m/2)/ y^) are dimensionless reciprocal lattice vectors, and Qm,n = exp[— 7r(g(^ + gy)]/(3A- 
In particular, 



7(o,o|o,o) = ^e„,„ = i 

m,n 

J(0,fc|0,fc) =^e„,„exp 

V 



|7(0,0|fc,-fc)| 



477^3 
^ 0m,n exp 



k"^ + 0.0029fc^ + 



ky ikx ) {igx 9v) 



1 - 



277^3 



r + o.oo32r + 



k^O. 



(D2) 

(D3) 

(D4) 
(D5) 

(D6) 



Within the BZ, it can be shown that < e_ < 0.59 < e+ < 1. The asymptotic behavior of e± can be shown to be 
at small k: 



e+{k) « 2 
e-ik) 



V 



k^ + O.OOQQr + • • 



^66/3^fc4 = 0.0026fc* + . 



(27r)2 



(D7) 
(D8) 



It is instructive to plot the contour of the function e±{k) for a triangular lattice configuration to illustrate its 
symmetries and the first BZ (see Fig. D). The maxima and minima are denoted by the light and dark shades 
respectively. The reciprocal lattice vector(RLV) points are marked either by the maxima of e+ or the minima of 
e_. The equal-sided hexagonal BZ can be constructed by joining together the six minima (maxima) surrounding the 
central RLV in the function e+ (e_). It is easy to see that approximating the two dimensional integration by a circular 
BZ underestimates the integrand involving e_ {k) as the function has spikes at the corners of the Brillouin zone. 
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FIG. 9. The contour plot of the functions (a)e+(fc) and (b)e-(fe) for the triangular configuration. 
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